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The Veneziano ghost field has been proposed as an alternative source of dark energy whose energy 
density is consistent with the cosmological observations. In this model, the energy density of QCD 
ghost field is expressed in terms of QCD degrees of freedom at zero temperature. We extend this 
model to finite temperature to search the model predictions from late time to early universe. We 
depict the variations of QCD parameters entering the calculations, dark energy density, equation 
of state, Hubble and deceleration parameters on temperature from zero to a critical temperature. 

We compare our results with the observations and theoretical predictions existing at different eras. 

It is found that this model safely defines the universe from quark condensation up to now and its 
predictions are not in tension with those of the standard cosmology. The EoS parameter of dark 
energy is dynamical and evolves from —1/3 in the presence of radiation to —1 at late time. The 
finite temperature ghost dark energy predictions on the Hubble parameter well fit to those of ACDM 
and observations at late time. 

PACS numbers: 95.36,+x, 98.80.Es, 12.38.-t, ll.10.Wx, 98.80.-k 


I. INTRODUCTION 

The universe is expanding at an increasing rate supported by various observations such as Supernova type 1A 
explosions, cosmic microwave background (CMB) radiation and baryon acoustic oscillations (BAO) jTHs]. There is a 
need for a kind of energy to fill roughly 75% of the universe that causes the late time accelerated expansion. The 
A-cold dark matter (CDM) model is currently the best cosmological model explaining this expansion. It is assumed 
that the cosmological constant, A, may arise from vacuum fluctuations. However, there is a large (at least sixty orders 
of magnitude) discrepancy between the predicted energy density of the vacuum in particle physics (of the order M 4 
with M p = fts 10 19 GeV being the Planck mass) and the energy density of the cosmological constant obtained 

from fitting the ACDM model predictions to observations, i.e. p odserved = (2.3 x 10 -3 eV) 4 [U, 3^1 ■ This is called the 
cosmological constant problem and many models have been proposed to overcome this problem. One of those that 
relates the vacuum energy to the QCD vacuum has been proposed by F. R. Urban and A. R. Zhitnitsky EMI at which 
Veneziano ghost field was firstly considered as a candidate for the late time acceleration. Veneziano first proposed 
this field as a ghost by putting a minus sign into the propagator with the aim of solving the U(1)a problem. This field 
is called ghost, since it describes the long range interactions of QCD [uj (for a review see Ref. |l4j)- In the model 
proposed by Urban and Zhitnitsky, the QCD vacuum energy is related to the fundamental QCD parameters and it is 
in the order of p vac ~ Hm q (qq) /m^' ~ (4.3 x 10” 3 eV ) 4 with H,m q , (qq) and m v ' being the Hubble parameter, light 
quark mass, light quark condensate and mass of the rj' meson, respectively. This energy density is in the same order 
of magnitude with the observations implying that the vacuum energy can be considered as QCD vacuum. The idea 
of using QCD in cosmology is not new; such that ample knowledge of QCD at finite baryon density and temperature 
has been used to understand a wide range of phenomena in cosmology such as the structure of neutron stars, Big 
Bang Nucleosynthesis (BBN) and so on. In Ref. 15], it is shown that different phases of QCD at finite temperature 
and density lead to interesting effects. In Ref. lfj, it is stated that the gravity may be a low energy effective theory 
of QCD, instead of being a fundamental interaction. It is believed that at very high temperatures (higher than the 
critical temperature), the quark condensates are not formed and the universe consists of quark-gluon plasma (QGP). 
In Ref. E3 the equation of state (EoS) parameter of QGP is investigated. 

The standard calculations in QCD are generally performed in Minkowski space-time with topological susceptibility 
X = 0 for zero scalar curvature. In Refs. EMI,' to investigate the QCD nature of dark energy, it is realized that if 
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our universe is embedded on a nontrivial finite manifold, such as a torus, the energy density should be proportional 
to the deviation from the Minkowski space-time. In this model, the topological susceptibility x is defined as 

X = iJ d 4 xe iqx (0 | T{Q{x),Q(0)} | 0), (1) 

where Q = is the topological charge density and 

K ' J = (d x A% + | f abc A\A%) , ( 2 ) 

with A“ being the conventional QCD gluon fields and g = y/4na a . Although Q is a total derivative, the standard 
Witten-Veneziano solution of U(1)a problem implies that x does not vanish. This causes an unphysical pole at zero 
momentum in the correlation function of K^. Using (0 | | ghost ) = Awith A ym and being the decay 

constant and the polarization vector, respectively; as well as e^e,, = —one can write 

lim i f d A xe iqx (0 | T{Q(x),Q( 0)} | 0) = -^-A YMq^q v , (3) 

9^0 J <r 

where is the propagator of ghost field and the minus sign is considered to solve the U(1)a problem By 

introducing a single light quark with mass m q and the matrix element of the rf field via (0 | K^ | rf) = =q^ one 
gets 

X = lim i [ d 4 xe^(0 | T{Q(x), Q(0)} | 0) = ^7°^ , (4) 

J (q z - m‘,) 

where m 2 , = + X 2 ,/N c is the mass squared of the physical rf field. On the other hand, the relevant Ward identity 

for QCD with light quarks 


X = i y d x e iqx {0 | T{Q(x), Q(0)} | 0) = m g (qq) + 0{m 2 q ), 


(5) 


is satisfied. By using these relations and m^f 2 = —4 m q (qq), the famous Witten-Veneziano relation, 4A^ M = f 2 m 2 ,, 
is obtained. The deviation mentioned above is related to the finite size (L) of the manifold and AA ym = = H. 

The topological susceptibility is related to the energy density through [(J jl8| 


X = 


d 2 e V ac{0) 

80 


e=o, 


( 6 ) 


where 6 is the angle of QCD. The corrections due to very large but finite size of the manifold are small. It is important 
to mention that if one knows the 0 dependence of vacuum energy, one can compute the energy mismatch that arises 
in theory between the infinite Minkowski and finite compact space-times. This effect is entirely due to the ghost and 
only much smaller corrections arise from all other QCD fields. From fl])-([n]), one can write (for details, see Ref. 1) 


d 2 e va c(6) 

80 2 


10=0 


= -Ax = A 


m 0 *YM 


~ —c- 


l r)i 


b 0 /X YM 


v 


2 N f H , , 

-c- J — | m q (qq) < 0. 

m v r 


(7) 


The 0 dependence of vacuum energy at 0 <C 1 and for Nf quark flavors with equal mass is given as e vac = | 

m g (qq) \ cos in Ref. ((I, HH,[2(| , such that d 2 e vac (0) = Note that all contributions from the gluon condensates 

and condensates from the Higgs field and etc., cancel out in the subtraction as they appear with equal magnitude. 
From (0 one can write 


Pg 


. 2 N f H . ... 

= Ae„ ac = c- | m q (qq) \ . 


l r\! 


( 8 ) 


This energy density in an expanding background, the Friedmann-Robertson-Walker (FRW) space-time, is analyzed 
in Ref. [l2j]. The global topology could be a torus and FRW metric still locally describes the space-time. The universe 
may have a non-trivial topology and there are different searches on this, such as the matched circle test. It must be 
mentioned that these searches yield no detection of a compact topology pll - [23| . 
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In Ref. jlH, Urban and Zhitnitsky use the QCD vacuum energy in FRW space-time to calculate some cosmological 
parameters and compare the model predictions with those of the ACDM and cosmological observations. In the 
calculations, they consider the QCD parameters at zero temperature (late time) to investigate the evolutions of EoS 
parameter of dark energy and Hubble parameter at low redshifts. In the present study, we extend those calculations by 
considering the QCD parameters at finite temperature. By increasing the temperature we can go from late time to the 
early universe and look for the variations of QCD vacuum, and as a result, variations of the cosmological parameters 
with respect to time (temperature). In particular, we use the temperature dependent m^ (T) and (qq)(T) to investigate 
the ghost energy density parameter (f l g ), EoS parameter of ghost dark energy (ui g ), total EoS parameter (uJtot), Hubble 
(H), and deceleration ( q ) parameters as a function of e-folding N and redshift z. The results are compared with those 
of zero temperature existing in the literature as well as the ACDM. 

The outline of the paper is as follows. We present some details of the model at finite temperature and its modifica¬ 
tions on the energy conservation in the following Section. In Section III, we discuss the cosmological parameters and 
their time evolution. The last section is devoted to the concluding remarks. 


II. THE GHOST ENERGY DENSITY AT FINITE TEMPERATURE 

As previously discussed, the energy density of QCD ghost field, p g , can be related to the Hubble parameter at late 
time as 


p g = aoH 0 , 


with 


ao 


2cN F \m q (qq) 0 \ 

m v ' 


(9) 


( 10 ) 


where (qq)o = (—240 x 10 6 eV) 3 is the light quark condensate at zero temperature, m q = 3.5 x 10 3 MeV is the 
average light quark mass (up and down), m v ' = 957.78 ± 0.06 MeV f24| is the mass of rf meson at vacuum, c is the 
speed of light, Hq is today’s value of the Hubble parameter (H = ^) and Np denotes the number of flavors. We use 
the natural units 87 tG = jp = c = 1. 

The parameter a o is constant at zero temperature. By setting Np = 2 and putting all other values, the energy 
density of QCD ghost field is found as p g ~ (4.3 x 10 -3 eU ) 4 which leads to a consistent prediction with the 
observations at late time. To investigate the dependence of a on temperature, we need to know the dependence of the 
quark condensate {qq) and on temperature. In the present study, we use the parametrization for the behavior of 
the quark condensate in terms of the temperature given in Ref. |25j which can be fitted to the following function: 


iw)(T) = {qq) o 


( _ 1 _ 

l l _|_ e 18.10042(1.84692T 2 +4.99216T-1) 


( GeV 3 ). 


( 11 ) 


This parametrization reproduces the lattice QCD predictions on the temperature-dependent light quark condensates 

EEE3- 

The behavior of rf meson mass with respect to temperature is studied in Ref. [28j and can be parametrized in 
terms of temperature as 


m V '{T) = 0.958 - 0.082T + 6.127T 2 - 79.287 T 3 {GeV), (12) 

where we took the critical temperature T c = 197 MeV (see also Ref. @)- Using the above parametrizations, we plot 
the dependence of the quark condensate, my(T'), and a{T) in terms of temperature up to T = 220 MeV in figure 1. 
From this figure we see that {qq) and m^ are constant up to some temperatures after which they start to diminish, 
drastically. The parameter a also remains unchanged up to roughly 0.16 GeV , after which it immediately falls to 
zero. We shall remark that the a parameter is barely constant for the late time as seen from the figure. To simplify 
the calculations, we use the e-folding (redshift) instead of temperature in the following. The relation between the 
scale factor, redshift and e-folding is 


Ta 0 e 

v — = — = l + z = - 
i o a e 


N 0 

AT 


(13) 


where N = In a and subindex “0” stands for today’s values. The parameter p takes the values p = (x) 3 f° r before 
the electron formation (T > m e ) and p = 1 for the first electron formation up to now (T < m e ). The dependence 
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FIG. 1: The dependence of (qq }, m v i and a on temperature. 


z (redshift) z (redshift) z (redshift) 



N (e-folding) N (e-folding) N(e-folding) 


FIG. 2: The (qq), m v i and a as a function of N (e-folding)/z (redshift). 


of the parameters under consideration on N/z is shown in figure 2. From this figure it is clear that the a remains 
unchanged up to N = —26.5, after which it starts to grow then immediately falls near to N = —28. 

Before delving into details, we would like to overview the constant QCD ghost implications in cosmology and 
compare it with the temperature dependent QCD vacuum. The Friedmann equation 

3H 2 Mp = a 0 H , (14) 

gives the exact solution of the scale factor for constant a as 

a(t)=a 0 e™t t . (15) 

The constant a at late time also implies w g = — 1 for the ghost dark energy and this gives the de-Sitter type expansion. 
When we take the temperature into account, the solution of Friedmann equation becomes more complicated. The 
continuity equation 


Pg +3(1 +Wg)Hp g = 0, 


(16) 


is naturally modified to 

(aH + aH) + 3(1 + w g )Hp g = 0, (17) 

at finite temperature. The extra aH term may be considered as an interaction term between the ghost field and other 
components of the universe. The interaction type is determined by the temperature dependence of QCD parameters 
but this is not investigated in the present paper. The cosmic evolution of the universe in chronological order from Big 
Bang up to now is inflation, radiation dominated era, matter dominated era, recombination, structure formation and 
late time accelerated expansion. In this paper, from now till the radiation-matter equality (T « 3 x 10 4 K, t ss 10 4 
yr, z rm = 2740, N rm = —9.3), the QCD ghost field and matter are considered to fill the universe and the radiation is 
safely neglected. After the radiation-matter equality, the radiation dominates over matter and the QCD vacuum and 
radiation are considered to fill the universe. 
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III. ANALYSIS OF SOME COSMOLOGICAL PARAMETERS AT FINITE TEMPERATURE 


A. Energy density parameter of ghost dark energy 


At late time (the ghost field and matter dominated era), the Friedmann equation is written as 

3 H 2 M 2 p =p g +p m , (18) 

which satisfies the continuity equations 

Pg + 3 (1 + Wg) Hp g = 0, and p m + 3Hp m = 0. (19) 

Dividing both sides of (11811 by the critical energy density p c = 3 H 2 Mp, we obtain f l g + fl m = 1 where is the 

dimensionless energy density parameter of x. The continuity equation for matter can be written in terms of f l m as 


and in terms of fl ff as 


fim + ( 2— + 3 H j fl m — 0, 


-tig + ( 2— + 3 H ) (1 — fl g ) — 0. 


The rate of change in Hubble parameter is related to fl g as 

H a tl g 
H = a ~ fy 

hence, m becomes 

ti g (n g - 2 ) + f 2 ^ + 3 hj (i - n g )n g = o. 


Let us rewrite (E3l) in terms of e-folding as 


Og(Og 


2 ) + 



(1 Hg)Og — 0, 


( 20 ) 

( 21 ) 

( 22 ) 

(23) 

(24) 


where ' denotes the derivative with respect to e-folding. 

For radiation and ghost field dominated era, considering the continuity equation for radiation, p r + AHp r = 0, from 
a similar manner, we obtain 

O;(Hg-2) + (2^+4)(l-Q 9 )f2g = 0. (25) 

The fig versus N for the two eras, obtained by the numerical solving of (l24l) and (l25l) , is depicted in figure [3] The 
vertical line on this graph shows the radiation-matter equality. From figure [3l we see that the ghost dark energy 
density at T yf 0 evolves exactly the same as the zero temperature case (a = const.). From this figure, we also see 
that there is no role of the ghost dark energy for N < —3 and it is produced when N > —3 in the universe. In Ref. 
[30j , the fig is analyzed in the interval —4 < N < 2 at zero temperature by considering no interaction between the 
ghost dark energy and CDM. We extend the interval up to N ~ —30. Our result (a ^ const.) is in a good agreement 
with the predictions of Ref. [30j] for fl 9 versus N in the interval —4 < N < 2. 


B. The EoS parameter of ghost dark energy 

For late time, rewriting the continuity equation for ghost dark energy csd in terms of fl g gives 

fig + | 2— + 3H(1 + Wg) J fig = 0. 


(26) 
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z (redshift) 



FIG. 3: The dependence of ghost dark energy density parameter on N/z. 


Using this equation together with (l2ll) and (l22l) . we get the EoS parameter as 


w 9, i = 


(2 - n g ) 


-1 - 


2 a' 

3 a 


(27) 


where the subindex 1 refers to the late time. Before the radiation-matter equality, in the presence of radiation and 
ghost dark energy, the EoS parameter reads 


W 9 ,2 


1 

3(2 — Ug) 




(28) 


with the subindex 2 being representing the radiation and ghost dark energy dominated era. We plot the dependence 
of w g on N/z for the two eras in the left panel of figure [4] We also would like to calculate the EoS parameter of the 
total fluid for two above mentioned eras. Since w m = 0, from the radiation-matter equality up to the finite future, 
N = 2, we have 


'Mtot - ^(/,lUg. (29) 

From QCD phase transition ( N qc d ~ —27) up to radiation-matter equality ( N rm = —9.3), we can also write 


Wtot - 'aJg^Ug + (1 Ug)lC r , (30) 

where w r = 1/3. We plot the dependence of w to t on N/z for the two eras in the right panel of figured] 

Before the description of the results presented in figure |4] and discussion on the evolution of EoS parameter with 
respect to e-folding, we shall remind that the ACDM defines a fiat universe with a cosmological constant (w\ = — 1 ) 
and wCDM extends this model to allow the EoS parameter to be different than —1. The Planck collaboration gives 
w = — 1.13lo!25 from combined Planck+WP+highL+BAO data (3l| and BOSS collaboration gives w = —0.97 ± 0.05 
from the most recent combined Planck+BAO+CMB data (32). In figured! when going from early to late time, we 
see that our model predicts w g = —1 near to the critical temperature but it immediately increases to w g = — ^ 
and remains unchanged up to the radiation-matter equality for a ^ const. This means that the ghost dark energy 
behaves as a cosmic string in the presence of radiation. However, in this interval, the QCD ghost is dominated by the 
radiation (see the right panel) and thus we can safely say that the QCD gho st does not change the BBN predictions. 
Our results on the w g and wtot are consistent with predictions of Ref. |30| for non-interacting case existing in the 
interval —4 < N <2. 

After radiation matter equality, the ghost field EoS parameter becomes roughly — i in the presence of matter and 
remains unchanged up to IV ~ —1, after which it evolves to —1 for now. In this era, the total EoS parameter is 
zero referring to the deceleration in expansion up to N ~ —3 then it starts to drop to —1 representing the late time 
accelerated expansion. The transition from the deceleration (Wtot > — |) to acceleration (Wtot < — 5 ) in expansion 
occurs at N = —0.51 (z ~ 0.6). This is a bit late compared to the ACDM prediction at which this transition takes 
place at N = —0.55 (z ~ 0.7). 

It is interesting that such variation of EoS parameter for dark energy at different eras has been predicted in some 
independent studies (33l-f36|. They have found w = — | for early time showing that forms of matter such as domain 
walls do exist, and it evolves to w = — 1 for late time. Similarly, in Ref. [36|, the authors re-parameterize the dark 
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FIG. 4: The EoS parameter w g versus N/z (left panel). The total EoS parameter wtot versus N/z (right panel). 


energy source and find that the dark energy source with a dynamical EoS parameter — | at early periods of universe 
and —1 today, matches slightly better than ACDM model to the recent observations. In this work, we find that the 
QCD ghost dark energy with a dynamical EoS parameter starts from the value — 4 at earlier times and goes to —1, 
behaving similar to the cosmological constant as time evolves. 


C. Hubble and deceleration parameters 

In ACDM model, the Hubble parameter is identified in terms of e-folding as 

—2 = Da,o + ^m,oe 3N + f^r,oe 4JV , (31) 

u o 

in FRW space-time for flat space-like sections. Here, the density parameters satisfy O g + O m + Q r = 1; H 0 = 
70.6 km/sec/Mpc , Da = 0.72, D m = 0.28 {H2| and Q r ,o = 2.47 x 10 -5 are considered. The Hubble parameter in terms 
of Wtot can be written as 

H + ^H 2 {1 + Wtot) = 0. (32) 

The variation of Hubble parameter (left panel) and its zoomed version (right panel) in terms of N/z are depicted in 


z (redshift) 
1.72 0 



3.48 


400 


z (redshift) 

1.72 0.65 


Ghost DE (ff=const.) 
ACDM 

Ghost DE (n^const.) 



N (e-folding) 


-1.0 -0.5 

N (e-folding) 


FIG. 5: Hubble_ parameter versus N/z (left panel). Zoomed version of H versus N/z together with data from the BOSS 
Collaboration ^3, [38] (right panel). 

figure [5j For comparison, in the same figure, we also show the variation of Hubble parameter in ACDM model as well 
as the variation of this parameter obtained at constant a by Ref. 0 , i-e- 

£ = n g>0 e-^+^ N » N + n m , 0 e~ 3N + fl r oe~ 4N . 

Ho 


( 33 ) 



















We also add the most recent data H3, [Hj to figure [5] 

From this figure, we see that our model’s prediction (a ^ const.) on variation of H is very close to that of the 
ACDM model and very well fits to the observational data. The behavior of H versus N for a = const, obtained in 
Ref. pjJ given in the interval —1.2 < N < 0 considerably differs from the others in the interval —1.2 < N < —0.2, 
although it predicts the same H with other models and observational data at N = 0. From this figure, we also see 
that, when going from the late to early time, the Hubble parameter drastically increases in terms of e-folding. 

Another cosmological parameter is the deceleration parameter, defined as 


. H 1 + 3w to t 
2 


(34) 


The deceleration parameter q versus N/z for late time is depicted in figure [6] From this figure, we read that the 


z (redshift) 

2.32 1.7 1.22 0.82 0.5 0.2 0 



FIG. 6: Deceleration parameter versus N/z. 

behavior of q with respect to N at late time shows roughly good consistency with the predictions of ACDM. It is 
well known that the transition from positive to negative values of q refers to the transition from deceleration to the 
present acceleration in expansion. As seen from figure EH in the ACDM model, this transition occurs at N ~ —0.55 
(z ~ 0.7), while in our model, the acceleration starts a bit later at nearly N ~ —0.51(z ~ 0.6). 


IV. CONCLUSIONS 

The ghost dark energy model proposed by Urban and Zhitnitsky @43 at zero temperature gives the energy density 
compatible with the cosmological observations at late time. We extended their model from late to early periods of 
the universe and searched for the variations of some cosmological observables in terms of temperature by considering 
the temperature-dependent QCD parameters. We discussed the variations of the energy density parameter of ghost 
dark energy, the EoS parameters of the ghost and total fluid, the Hubble and deceleration parameters with respect 
to N/z. 

It has been found that the ghost dark energy plays no role at early universe and it is produced after N > —3. 
The lack of ghost dark energy before N > —3 is in agreement with the predictions of the standard Big Bang on the 
decelerated expansionary phases (matter and radiation dominated eras). Our prediction on S2 s is in a good consistency 
with the prediction of Ref. [3 existing in the interval —4 < N < 2 for constant a and non-interacting case. 

The model predicts the EoS parameter of ghost dark energy to be w 9 = —1 near to the critical temperature, 
although it drastically increases to w g = — 4 referring to the behavior as a cosmic string in the presence of radiation. 
It remains unchanged up to the radiation matter equality. In this interval, the total EoS parameter of the fluid shows 
that the QCD ghost is dominated by the radiation and the predictions stay consistent with the BBN’s and so with the 
values for the abundance of light elements. After radiation matter equality, the ghost field EoS parameter becomes 
roughly —4 in the presence of matter and remains unchanged up to N ~ — 1, after which it evolves to —1 for late time. 
After radiation-matter equality, the total EoS parameter becomes zero representing the deceleration in expansion up 
to N ~ —3 , after which it starts to drop to —1 referring to late time accelerated expansion. The transition from the 
deceleration to acceleration in expansion happens a bit later in our model compared to the ACDM. The late time 
predictions for w g and Wtot. are in good agreement with those of Ref. [30j| for constant a and non-interacting case. 

In the case of Hubble parameter, our prediction is consistent with that of the ACDM and better fits to the 
recent observational data at late time compared to the prediction of Ref. [3 f° r constant a existing in the interval 
— 1.2 < N < —0.2. All models have the same predictions at N = 0. 
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For the deceleration parameter, our model predicts that the change in sign of this parameter (transition from the 
decelerated expansionary to the accelerated expansionary phase) occurs a bit later than that of the ACDM model 
although these models have roughly the same predictions on the behavior of q with respect to N at late time. 

The obtained results at finite temperature in the present work point out that the QCD ghost vacuum can be still 
a valid candidate to dark energy. Interestingly, this dark energy source has a dynamical equation of state parameter 
equal to —1/3 at the early universe and —1 today, behaving similar to the cosmological constant at late time. This 
dynamical property of EoS parameter can be checked by fitting to observations. 
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